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Abstract 

The material tensor of linear response in electrodynamics is constructed 
out of products of two symmetric second rank tensor fields which in the 
approximation of geometrical optics and for uniaxial symmetry reduce 
to "optical" metrics, describing the phenomenon of birefringence. This 
representation is interpreted in the context of an underlying internal ge- 
ometrical structure according to which the symmetric tensor fields are 
vectorial elements of an associated two-dimensional space. 
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1 Introduction 

Electrodynamics in media is a broadly investigated subject with a variety of 
different facets As far as the propagation of waves is concerned, the 

approximation of geometrical optics has played a major role. Within this ap- 
proach Maxwell's equations are transformed into a set of algebraic equations 
for the wave vector. In a relativistic context (see 01 EJ El El ) j a specific feature 
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of the wave propagation in media is the appearance of an "optical" metric, in- 
troduced by Gordon With respect to this quantity light in isotropic media 
propagates as in vacuum. Light rays follow geodesies [5] and the wave vector is 
a null vector of this metric. The optical metric is characterized by the refrac- 
tion index of the medium which is composed of the coefficients of dielectricity 
and magnetic permeability. The latter quantities, in turn, enter the material 
tensor that relates the excitation (induction) tensor to the field strength tensor. 
For the motivation of this paper it is essential that this happens in a manner 
which allows one to represent the fourth rank material tensor as a combination 
of products of the second rank optical metric. Written in terms of the optical 
metric, the material tensor for an isotropic medium has the same structure as 
the "material" tensor for the vacuum, only that the spacetime metric within the 
latter is replaced by the optical metric [S] . Here we ask to what extent this type 
of representation can be generalized to anisotropic media. This is anything but 
straightforward since light propagation in anisotropic media is accompanied by 
the phenomenon of birefringence. A hint here is the fact that in the simplest 
anisotropic case, the case of uniaxial symmetry, birefringence can be described 
by two Lorentzian optical metrics [7| (under more general circumstances they 
have to be replaced by Finsler metrics). Given the mentioned representation of 
the material tensor in the isotropic case, this raises the question whether there 
exists a bi-metric representation of this tensor in anisotropic media. While the 
problems of birefringence and bi-metricity have been considered for non-linear 
electrodynamics (birefringence as a non-linear phenomenon) |51 110|. we are not 
aware of an exlicit representation of the material tensor in terms of two optical 
metrics in the linear case, describing the so-called ordinary and extraordinary 
waves in uniaxial media. We shall demonstrate that the material tensor can 
indeed be written as a combination of products of those two optical metrics 
which characterize the light cone structure for uniaxial symmetry. A construc- 
tion of the material tensor out of two symmetric second rank tensor fields is 
even possible in general anisotropic media, but the interpretation of these quan- 
tities as optical metrics is restricted to the uniaxial case. We also show that the 
transition between different representations of the material tensor in terms of 
different pairs of second rank tensor fields is governed by invariance properties 
in an associated two-dimensional internal vector space. 

The paper is organized as follows. In section|3we consider Maxwell's equa- 
tions and the general constitutive laws which mediate between the tensors of 
excitation (induction) and field strength. We decompose the material tensor 
with respect to contribution parallel and orthogonal to the four velocity of the 
medium. We recall the role of an optical metric in isotropic media and provide 
different representations of the material tensor in terms of this metric. Section 
is a collection of basic relations of geometrical optics. Section which is the 
main section of the paper, is devoted to the construction of two symmetric sec- 
ond rank tensor field in terms of which a representation of the material tensor 
is obtained. Novel geometric features of this representation are pointed out. In 
the approximation of geometrical optics Maxwell's equations reduce to a fourth 
order algebraic equation for the wavevector, the "extended" Fresnel equation 
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(cf. For the case of uniaxial symmetry this equation is shown to factorize 

where the two symmetric tensor fields reduce to optical metrics for ordinary and 
extraordinary wave propagation. A summary of the paper is given in section |3J 



2 Maxwell's equations and constitutive laws 
2.1 General theory 

The covariant, source free Maxwell equations in continuous electromagnetic me- 
dia are 

V fc F * lk = and V k H lk = , (1) 

where 

F* lk = ^=e lkLs Fi s (2) 

is the dual to the field strength tensor F lk and H lk is the induction (or ex- 
citation) tensor. Both F lk and H lk are antisymmetric, i.e., F lk = —F kl and 
H — —H kl . Here, Vfc denotes the covariant derivative with respect to the 
spacetime metric (with determinant g and signature H ) and — =e lMs 

is the Levi-Civita (pseudo-) tensor where e tkls is completely antisymmetric with 

,0123 _ , _ -i 
e — —6oi23 — J-' 

The set of equations (JTJ has to be completed by constitutive relations, link- 
ing the electromagnetic inductions to the field strengths. Restricting ourselves 
to linear dielectric, permeable media, the most general structure for the corre- 
sponding constitutive laws is 

H ik = c ik ™F mn , (3) 

where C mn is the material tensor which characterizes the electromagnetic 
properties of the medium in the linear response approximation. This tensor has 
the symmetries 



Consequently, it has 21 independent components. (See, however, [1 1 j for an 
approach without the symmetry C tkmn = C mmk ). Assuming that the medium 
is characterized by a four- velocity field U l , normalized by U l Ui = 1, we may 
define the vectors D\ H\E\ and B l as Ql] 

D i = H ik Uk ^ H i = H *ik Uk ( E i = F ik Uk ; B t = F *ik Uk _ ( 5 ) 
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These vectors are orthogonal to the four- velocity vector U l , 

D l Ui = = E% , H l Ui = = B% . (6) 
The field strength tensor can be decomposed according to 

F mn = SPJ n E p U q ~ r, mnl B l , (7) 
where (5^„ is the generalized 4-indices 5— Kronecker tensor 

cik ci rk efe 

mn m n ^n^m 

and rj mn i is an antisymmetric tensor orthogonal to U 1 , defined as 

Vmnl = ^e mnls U s , ri m = -L=e lkls U s . (8) 

A useful identity is 

- rf^r, mnp = 6% l ns U t U s = A l m A k n - A> n A k m , (9) 

where 

gikl = gi gkgi + S k gi gi + gi gi § k _ gi gi § k _ gi gkgi _ gk gi gi 

and 

A lfc = g lk - U l U k . (10) 
Upon contraction, equation J^J yields the further identity 

\v M m m = -sLuiU' = -A l m . (ii) 

2.2 Decomposition of the material tensor C tkmn 
We can uniquely decompose the tensor C lkmn as 

^jikmn jj_ ^imjj-kjj-n ^injjkjjin _|_ ^knjjijjm ^kmjjijjnl^ 

-I [rf k \U m Vl n - C/V; ro ) +T} lmn (U i v l k - 17*1/,*)] . (12) 
Here e™, {(J-~ 1 ) pq and v p m are defined as 

e im = 2C lkmn U k U n , 
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/ — 1 \ _ siikmn 

VP 1 Jp<? 'Ipik^ 'Imnq j 



VpikC ikmn U n . (13) 



The tensors Eik and (fi~ 1 )ik are symmetric, but v k , is in general non-symmetric. 
The dot denotes the position of the second index when lowered. These three 
tensors are spacelike, i.e., they are orthogonal to U l , 

s ik U k = 0, (fi-% k U k = , v k \J l = = ^ k U k • (14) 

The total trace of the C %kmn tensor is equal to 

C nkmn g im g k n = s k + {^t . (15) 

With the help of the definitions Q , © and the decomposition l|12(l , we obtain 
the linear laws 

D l = e m E m - B l v ] i . and H t = {^% m B m + u t m E m . (16) 

Obviously, the tensors e lm and pt pq are the four-dimensional analogues of the 
dielectricity tensor and the magnetic permeability tensor, respectively. The 
tensor i/ p m describes magneto-electric cross effects. Thus, the 21 components of 
Qikmn cons i s t f the 6 components of e lm , the 6 components of (fJ*~ 1 ) pq and the 9 
components of v™.. For media without cross effects the number of independent 
components reduces to 12. 

2.3 Special cases 

2.3.1 Vacuum 

The vacuum is characterized by 

ti = -t - g g t mn , (17) 
which corresponds to a "material" tensor 

(vac) - 2 \ g 9 ~ 9 9 ) ■ l-LSJ 

Under this condition the quantities l|13|l reduce to 

e* k = (^ 1 y k = A lk , v\ k = . (19) 

Alternatively, introducing the latter relations into Eq. Ijl2(l , we obtain the special 
case (fT£|) . 

2.3.2 Isotropic medium 

A spatially isotropic medium is characterized by 

D l = eE l and H t = -B, , (20) 
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with two phenomenological scalars, the dielectricity e and the permeability \i. 
The relations l|2U[l are the special case of the general structure for 

e lk = eA lk , (^ 1 )' lk = -A lk , v k . = . (21) 

M 

For e and fi equal to unity, we recover the vacuum case Ijl9|l . Introducing the 
expressions (J2U for e lk , (^ 1 ) tk and v k . into JT^) one obtains the material tensor 

ikmn _ ( Am „kn An .fem\ i ( '-M ^ \ / Am jjkjjn 



s-iikmn ( Am „kn An „km\ , I C P \ { AmjjkjT 

°(isotr) - 2~ [9 9 ~99 ) + I 2 ^ J ^ ^ t> 



-g m U K U m + g Kn U l U m - g Km U l U n ) . (22) 
A convenient representation is obtained in terms of the second rank tensor 

where n is the refractive index, defined by n 2 = efi. The inverse tensor is 

= n 2 A krn + U k U m , (24) 



* _ 2 
9km — n 



9fcm + I :3 - 1 ) C4C/„ 



such that 

9* lk 9* km = 51 . (25) 

Obviously, we have 

g» k U k = U\ g* ik UiU k = l, g* km U m = U k , g* km U k U m = l, (26) 

i.e., the tensors g* %k and g* km preserve the norm of the four- velocity vector. 
(Note, however, that generally indices are raised and lowered with the metric 
tensor gi k .) Furthermore, the relation 

An 

g* ln = ^ + U l U n (27) 
en z 

for the isotropic case e %k = eA %k is valid. With the help of both g %k and g* lk 
the tensor IL'l't can be written as 

77 f 77 

°(isotr) - 2^ Uf -9 ~5 5 J ~ ( n 2 _ x) U-9 ~5 J (3 ~9 ) 

I in *in\ l km *fcm\l 1 /oo\ 

- [g -9 ){g -g )\} , (28) 

which will be useful for later reference. The most compact and elegant way, 
however, is the symmetric vacuum type form (cf. |BJ) 

4 

(isotr) — 2^ \9 9 ~99 ) ■ \^) 



G 



In such a form the material tensor C 1 ^™" for the isotropic medium, up to the 
factor has the structure of (|T5|) . with the metric g lk replaced by g* lk . 

3 Geometrical optics 

3.1 Basic relations 

We assume the field strength tensor to be given in terms of a four potential A m 

by 

-^mn ^ mA n ^ nA m ■ (30) 

Geometrical optics is based on the leading-order approximation in the eikonal 
derivatives of 

Ai = a,e 1 *, k m = V m $, (31) 

where $ is the phase, a; is a slowly varying amplitude and k m is a wave four- 
vector. In this approximation (for complex quantities always the real part has 
to be taken) 

F"mn — i (k m A n /c nJ 4 m ) (32) 
and Maxwell's equations reduce to 

k\C (k m A n — k n A m ) = . (33) 

3.2 Vacuum 

In the vacuum case Eq. takes the form 

k^g^hA^-A^^hk^^Q. (34) 

Contraction with Ui yields 

(IW) (g ln kiA n ) = (U t A l ) (g lm hk m ) . (35) 

For the Landau gauge UiA 1 = and Uik 1 ^ it follows that 

g ln kiA n = . (36) 

While this is identical to the Lorentz gauge, it is not imposed as a separate 
condition in the present context but it is a consequence of the condition UiA 1 — 

0. With H3(j|) the field equation (|34|) admits a nontrivial solution for the potential 
A 1 if and only if 

g lm hk m = , (37) 

1. e., the wave four vector k % has to be a null vector with respect to the metric 
g tk . Moreover, the curves to which k a is a tangent vector are geodesies |13) . 
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3.3 Isotropic medium 

With the help of the optical metric (|23|) and with H29[) the Maxwell equations 
for an isotropic medium can be written in the compact form 




Proceeding as in the vacuum case, projection along Ui leads to 

(U^) (g* ln hA n ) = (U t A l ) (g^kkrn) . (39) 

For UiA 1 = (Landau gauge) and Uik 1 ^ we get 

g* ln hA n = , (40) 

which, because of UiA 1 — coincides with (|36l) and the Maxwell equations l|38[) 
simplify to 

g* m A n g* lm hk rn = . (41) 
This equation admits a nontrivial solution for the potential, if and only if 

g* lm hk m = , (42) 

i.e. if the wave four- vector k m is a null vector with respect to g* m , which justifies 
the name optical metric 0] for this quantity. It was proven by Pham Mau Quan 
|S] (see also ^J) that the corresponding light curves are null geodesies of the 
metric g* k . 

4 Anisotropic media 

4.1 Representation of the material tensor 

In the isotropic case it turned out to be possible to find a metric g* ik in terms of 
which the matter tensor has the quasi- vacuum representation l|29|) . Moreover, 
g* lk is an optical metric. With respect to g* lk the wave vector is a null vector. 
Both for the representation of the matter tensor and for the light propagation in 
geometrical optics the transition from vacuum to an isotropic medium consists 
in replacing the spacetime metric by the optical metric g* lk . It is natural to ask, 
to what extent these features may be generalized also to anisotropic media. A 
useful hint here is the circumstance that for anisotropic media of the simplest 
kind, the uniaxial crystal, there exist two optical metrics which account for the 
phenomenon of birefringence (see, e.g., [7]). This seems to suggest a bi-metric 
representation of the matter tensor in such media. The construction of C lkmn 
out of two different symmetric second rank tensor fields is the subject of the 
present section. While these tensor fields cannot be expected to be optical 
metrics in the general case, for an uniaxial symmetry, however, a decomposition 
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of C lkmn in terms of two optical metrics will indeed be possible. With respect 
to each of them light propagation proceeds as in vacuum in a similar sense 
as described by l|42[) for isotropic media. To the best of our knowledge, this 
problem has not been addressed in the literature so far. Our strategy will be as 
follows. In a first step we simply assume the existence of two symmetric tensor 
fields as construction elements of the material tensor. This choice is not unique 
and we have to single out a suitable pair of symmetric tensors. This will be done 
with the help of a geometric interpretation according to which these tensors can 
be regarded as vectors in a two-dimensional space. The next step will then be 
the explicit construction of these tensor fields and, finally, the demonstration 
that they reduce to two different optical metrics in the appropriate limit. 

This strategy implies that we consider media for which the tensor C lkmn in 
l|12l) is reducible to the form 

C lkmn = JL G («)(/3) (g vm{a) 9 kn(l3) ~ 9 m(a) g km{ ^ , (43) 

where g lm ( a ) and g lm ^ are a pair of (for the moment unknown) symmetric, 
non-degenerate tensor fields. The indices a and (3 take the values a and b and 
the sum is over all possible combinations. By construction, the expression l|43[) 
satisfies the symmetry conditions Q. For a given tensor C lkmn with generally 
21 independent components the representation may be regarded as a set 
of algebraic equations to determine the 10 components of g zm w , the 10 compo- 
nents of g lmi ^> and, additionally, one of the three quantities G( )(o)/Ai G (&)(&)//*; 
and G( a )( b )//t = G^^/fi (we assume G( a \rm to be symmetric). Two of these 
quantities remain arbitrary, the corresponding freedom will be used later on. 
The representation (|43[) is a direct generalization of (|18(l and l|29(l for the cases 
of vacuum and isotropic media, respectively. Since we require for the isotropic 
limit g im ^ = g im ^ = g* im the decomposition ^ with fj, = % to hold, we 
obtain 

E G C«>(0 = 1 ( 44 ) 

as an additional condition. The choice of the symmetric tensors g lm ( a ) and 
g lm i b ) j n the decomposition (|4*3|l is not unique. A first guess could be (cf. [TT] . 
(D.1.80)) 

/-likmn \ /'^,* m (l) „kn(2) „i«(l) „fc™(2) i rl im{2) „fen(l) „*«(2) „km(l)\ 

°(1,2) ~ 4^t \ / ' 

(45) 

which corresponds to 

G (l)(l) = G (2)(2) = 0, G(l)(2) = G (2)(l) = 2> dct ( G (l,2)) = (46) 
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where det(G(i j 2)) is the determinant G(i)(i)G(2)(2) — G(i)<2)Gt2)(i) of the ma- 
trix G( ct )( / 3) for this case. So far we will not use, however a specific choice but 
further consider the general representation (|43|l . As any symmetric tensor, the 
quantities g lk ^ can be decomposed with respect to their components parallel 
and orthogonal to the four- velocity U % , 

where 

B {a) =g iHa) UiU k , V p{a) = A P g lk(a) U kl S pq(a) = A p g tk{a) A q k . (48) 

In case U l is lightlike with respect to one or both tensors g lm ( a \ the correspond- 
ing coefficients B^ vanish. We shall assume, however, that, analogously to the 
vacuum and isotropic media cases, the four velocity U l is timelike with respect 
to both tensors g tm ( a \ Under this condition the factors B^ in front of U l U k 
in H47|) are different from zero and can be absorbed into the coefficients G( ce )( l g). 
The decomposition (|47|) then reduces to 

g lk(a) = yljjk + £>i(a)jj* + V k( a )jji + S ik(a) ^ ^ 

with quantities V p ^ and S pq ^ redefined accordingly. 

The magneto-electric cross terms i/ m p in are determined by the coefficients 
V p ^ via 

V ™ P = h E G (*m^k (v k ^g im ^ - g km ^) . (50) 

M (")(/3) 

Obviously, a sufficient condition for excluding magneto-electric cross effects is 
= 0. (Note that even in the presence of cross effects one of the V p ^ 
can be transformed to zero by assuming U 1 to be an eigenvector of the corre- 
sponding g lk ( a }). In the following we shall restrict ourselves to V p ^ — and use 

g ik(a) _ jjijjk + S ik{ a ) ^ U t g lm{a) = U m . (51) 

Applying the structure (|51|l to l|43|l. the spatial tensors e lm and (/J.~ 1 ) P q in l|13fl 
can be written as 

e im = ~ { [G iaKa) + G {am ] 5*«W + [G (am + G mh) ] S"^} , (52) 

and 
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(M 



-l\pq 



(a)(a)0 y 'O + Lr(b)(b)0 w 6 



(53) 



respectively. In a next step we make use of the remaining two degrees of freedom, 
mentioned in the discussion following equation l(43|l ■ From (|52|l it is obvious that 
there exists a special configuration, for which e %m depends only on one of the 
spatially projected (cf. 1)49(1 ) metrics, say on <S"™( a ). This configuration, which 
we denote by (a) = (A) and (b) ~ (B), is characterized by the two relations 



G(A)(A) +G( B )(A) - 1 



G 



(A)(5) 



G(B)(B) — 



(54) 



which obviously satisfy 144|) . Introducing the abbreviation 7 = Gr^m 



G 



(B)(A) 



-det(G( J 4 B)), where the latter quantity denotes the determinant 
G(A)(A)G(B)(B) - G( A ){B)G(B){A) of the matrix G(a)(b) for this configuration, 
the material tensor 1431) takes the form 



s~iikmn 



2u \ 



gi m (A)gkn(A) _ 5 m(A) 5 fem(A) 



fgim{A) _ gim(B)^gkn(A) 



-,kn(B)) 



(.9' 



in( B )\(gkm{A) 



km{B)\ 



)]}, (55) 



which generalizes the representation (|28|l for the isotropic case. The factor 7 
generalizes the (positive) expression n "_ 1 in (|28|l and has to be positive as well 
(cf. subsection 14. 21 below) . The quantities l|52l) and (|53|l specify to 



gim _ Gjim(A) 



(56) 



and 



respectively. 

It is evident from (|51|l that l|56l) implies 

1 r 



- 7 (s irn ( A )-$ tm ( B A {$ kn ( A )^$ kn ( B )\ ; (57) 



(58) 



gim(A) _ Jjijjri 



independent of the value of the parameter 7. This parameter contributes to the 
tensor (^ 1 ) pq only. In turn, (|58|l provides us with an explicit expression for 
terms of the tensor e m , which in the isotropic limit reduces to the 
optical metric l|27|l . 
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4.2 Geometric implications 

The change between the explicit representations H45|) to (|55|l of the material 
tensor corresponds to a transformation (7 > 0) 

g ik(l) = gl k(A) + ^ ^*(A) _ g ik(B)^ (5Q) 



and 



or, inversely, 



and 



g ik(2) = g ik(A) _ ^ fik{A) _ g ik(B) ^ f (6Q) 



ik(A) 
y 2 



l(g lk(1) +g lK2) ) (61) 



g ik(B) = gl k(A) + ^ g ik(2) _ g ik(l)^ (g2) 

Transformations between different sets of symmetric tensors in the representa- 
tion (|43fl can be regarded as transformations in a two-dimensional space with a 
metric G( Q )(/5), which is symmetric and non-degenerate, i.e., 

G {a)(0) = G{/3){a), G = det(G( a )( i a)) ^ 0. (63) 

This geometric interpretation requires a Kronecker delta S&l to exist together 

with a contravariant tensor GOOC"), satisfying G^^G^^ — fi9p\- Further- 
more, in such a space there exists an anti-symmetric Levi-Civita (pseudo-) tensor 



E (a)(/3) = V\G\ £(<*)(/?), e (a)(/3) = -£(/3)(a), £(1)(2) = 1- (64) 

The point here is that the decomposition (|43|l is invariant with respect to trans- 
formations in this associated space. In order to demonstrate this, we consider 
linear transformations 

9 ' -\ a )9 , G {a y W - t {a yt {f3 y& (am , (65) 

and require that the scalar product G ( a )(p)g mnt " a ^ 'g lk ^ remains invariant, i.e., 

(~i „im(a)'ln(0)' /~i n im(a) ln(@) 

G( a )>(pyg g l*(a)(p)9 9 ■ (66) 
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Contraction with UiU m yields 

E G(<*)>m-9 lnW = E G (am g ln ^. (67) 
(«)W (<*)(£) 

Identifying now (a)' and (/3)' with (A) and (£?), respectively, as well as (a) and 
with (1) and (2), respectively upon using (|51Jl and we reproduce the 
transformation l|tjl|) . Similarly, the invariance 

irn(a)' ln(fi)' 771 „im(a) ln(0) {ao\ 

E{ a )>(f})'9 y ' 9 y - E {a )(f))9 y '9 y ' (68) 

via 

E E {a y W ^' = E ^m^r (69) 
(")'(«' («)(£>) 

provides us with the second transformation (|62|) . As a consistency check, we 
also realize that 

G'=I^| 2 G = 4 7 G, (70) 

which is the correct relations between the determinants G' = det(G(A,B)) and 
G = det(G(i, 2 ))- Here, g in (fTTTfl stands symbolically for the set g lk ^\ g lk ( 2 \ 
and g' stands for the set g tk ( A \ g tk ( B ) . 

The relations (|59l) - (|70|l constitute an underlying geometric structure for 
the representation of the material tensor in terms of second rank symmetric 
tensor fields. The coefficients Gy/m play the role of a metric in the associated 
two-dimensional space, the tensors g lk ^ are vectorial objects in this space. 



4.3 Maxwell's equations 

With the general expression Maxwell's equations are 

C amn hk m A n = 

^ ^ E G («)(/3) (/ m(Q) 9 lnW ~ 9 m[a) 9 lm(f3) ) hk m A n = 0. (71) 

Contracting these equations with Ui similarly as for the previous isotropic case, 
we obtain 

U m k m £ G (Q)(/3)5 te ^)fcM„ - C/"^n E 5 im(/3) ^fc™ ■ (72) 

(«)(« (")(/9) 
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With the Landau-gauge U n A n = and U n k n ^ 0, the relation generalizing l|40|l 
is 

g\: s) hA n = , g[ n cS) = G {am g ln ^ . (73) 

(«)(/9) 

The structure of this relation suggests to identify g£ S \ with one of the symmet- 
ric tensors g ln ( a \ let's say g ln ( a \ out of which the material tensor is constructed. 
This is equivalent to 

5(cff) = 9 ln(a) = [ G (a)(a) + G (b){a)] + [ G (a)(b) + G (b)(6)] s'™^ 

G (a)(a) + G (b)(a) = 1) G( a ) (ft) + G(ft) (ft) = . (74) 

It is remarkable, that this choice again results in the previously derived set of 
coefficients, characterized by |J52J|. The identification gife^ = g ln ^ in (j?4"|) is an 
alternative way to single out this option which we again denote by (A) and (B). 
In the present context we have the additional relation 

g ln ^hA n = , (75) 

which has the structure of a generalized Lorentz gauge. We emphasize again 
(see the remarks following (|36l) ) that this is not imposed as a condition here but 
that it follows from the Landau gauge U n A n = 0. 

4.4 Constructing g kn ( A ^ and g kn ^ 
4.4.1 General relations 

While the tensor g lk ^ A "i is given in terms of e lk by i|58|) , it is less straightforward 
to obtain a corresponding expression for g lk i B ) , Following [7j we consider the 
eigenvectors Xhy X l , 2 y and of the tensor e tk . These eigenvectors have the 
properties 

9ikX\ a) X( b) = -5(a)(b), gikX\ a) U k = 0, (76) 
*(i)*(i) + X U X U + X U X k = ~9 lk + U l U\ (77) 

Vkl X (a) X (b) = e(a)(6)(c)X( c) , (78) 

where (a), (&).... = (1), (2), (3) are tetrad indices and the double index (c) de- 
notes a sum. Relation (|78H is a consequence of the fact that the scalar function 
ViklXf^XkjX 1 ,^ is totally antisymmetric with respect to indices (a), (6), (c) and 
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has to be proportional to the Levi-Civita symbol £( a )(6)(c)- The dielectric per- 
meability can be decomposed according to 

^ = - £ M X (a) X (a) , 4 = £(1) + £(2) + £(3) , (79) 
(a) 

where the terms £( a ) denote the (positive) eigenvalues, corresponding to the 
eigenvector X 1 ,*. . Since the tensor e lk is orthogonal to the four- velocity vector 
U l , the corresponding eigenvalue £( ) is equal to zero, and the velocity does not 
appear in this decomposition. 

The tetrad components &( Q ) and of the wave four-vector and of the 
electromagnetic potential four-vector, respectively, are 

fc (a) sMfo, A (6) =A^ a) (80) 

and we may write 

V = 11^11%) -k {a) X\ a) , #=-^,4). (81) 

The tensor fields g lm ( Q ) are describable in terms of U l and X? v according to 

(cf. (53) 

3 m(Q) = C/ 4 C/ m + sg (6) ^ a) X$ , (82) 

where 

5 («)M =5( ™ X M X M , (83) 

and where we again have neglected the magneto-electric cross effects. In the 
following we use the first and second equations (|T3|) to construct the tensors 

g im(A) and g im(B)^ 

4.4.2 The tensor g lm ^ 

While g lmi < A } is already known to be given by (|58|l . we briefly sketch how it 
is alternatively derived in the tetrad formalism. Together with (|82|l the first 
equation (|13fl becomes 

2C ikmn U k U n = -s {a) 5^X{ a) X% y (84) 
With A 4fe = -d^^X^XL it takes the form 

^ X (a) X (b) G («)(/3) 5 (a) (f,) = - £ {a)& {a)(h)X la) X (b)- ( 85 ) 
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Contraction with the tetrad vectors yields 

E <W>Sgj« = -^(*)S {aKb) ■ (86) 

(a)(0) 

With the coefficients (I54f) we obtain 

5y=-As (a) ^W, (87) 
such that the first of the tensors is 

g im{A) = jjijjm + j^im (gg) 

This reproduces (|58|l . For the choice fi = -|y and for e lk — eA tk we recover the 
optical metric 127JI for the isotropic medium. 

4.4.3 The tensor 

Applying the tetrad representation Q82JI to the second equation yields 

- lv P *kC* kmn Vmnq = -(jjL-% m X^X( b \ (89) 

where — (/i _1 )( a )(t>) ar e the tetrad components of the magnetic permeability 
tensor. With l|78|l it transforms into 

G (a)08) 5 (a) d)S (j3) (I e («X<0 W e (&)(«9C/) = 2 A(M _1 )(o)(6) ■ (90) 

Here e( a )( c )(e) is the Levi-Civita symbol in three-dimensional space, tetrad in- 
dices being regarded as Euclidean. The equations ((971)1 form a system of 6 
quadratic equations for the 6 components of S^j b \ Let us restrict ourselves to 
the case of a diagonal {n~ 1 )( a )(b) tensor, i.e., to 

(A* _1 )(a)(6) = —S(a)(b) , M(a) > 0- (91) 

M(o) 

This ansatz allows us to search for an explicit solution of the basic equations 
P3)l for which the 5 ( ( ^ ) ) (b) are diagonal, i.e., 

S$ (6) =0, if (a) ^(6). (92) 

Using the structures Ij54(l of Gt a \{p\ and 1|) of (/x -1 ) (a)(6)) we obtain from H90JI 
the system of equations 

< j(2)(2) < j(3)(3), < j(2)(2) /> ,c(3)(3) r ,. , /, -1\ r ,2_ _ _ A 
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c(l)(l)c(3)(3) , o(l)(l). ic( 3 )( 3 ),-,o ~-V\ ,»,2. . 

6 (B) 6 (B) + b (B) ^ £ (3)+ b (B) ^(l) + (l-7 JA*£(1)£(3)- 

<,(1)(1) <,(2)(2) ,<j(l)(l) f „(2)(2) r n r o 

6 (b) 6 (B) + 6 (B) M £ (2)+^( B ) M£(i) + il-7 jM£(l)£(2) 
After the substitution 

<,(<»)(&) _ _ - r(a)(6) , y(a)(b) 



7M(3) ' 



-, (94) 
(95) 

(96) 



where the first term on the right-hand side coincides with <S/^\ from 1)87(1 . the 
system (|93 |l - l(95(l takes the form 



F (2)(2) y( 3)(3) = ^ £(2)£(3)M(l)) 



A 2 



(97) 

(98) 
(99) 

(100) 

M(a)M£(l)£(2)£(3) 

In the isotropic limit with £( a ) = e, (j,^ = [i and (1 = l/(en 2 ), the M( Q ) 
quantities vanish identically. Solving the symmetric system (|97ll - <|99[) leads to 
the following solutions of the system l|9"3" |l -l|95 [l : 



y(D(i) r (3)(3) = ^ (1)£(3)M(2)i 
y (1)(1) y (2)(2) = ^ (1 ) £ (2)M (3) , 



where 



M, 



(a) 



1 



£(a) 



5, 



(i)(i) 



5, 



(2)(2) 



-M£(2) 



5, 



(B) - -M£(3) 



'(B) 

,(3)(3) 
(B) 

Thus, the second symmetric tensor is 



1± 



1± 



1± 



'M( 2) M (3) 
7% 



7 M (2) 



/M (1) M (2) 
7 M (3) 



M (1) M (2) M (3) 



(3) 

£ 

(»)=(!) 



Ha) 



(a) 



A (a) A (a) 



(101) 

(102) 
(103) 

(104) 



The first two terms on the right-hand side coincide with g lm ^ from (|88|) . In 
the isotropic limit we have g lm ( B ) = g lm ( A ) = g* lm . The parameter 7 remains 
undetermined. Using the explicit expressions (|88|l and (|104f> in i|55|) it is obvious 
that C lkmn is independent of 7. A proper choice of 7 will be essential, however, 
to make g lm ^ an optical metric in special cases to be discussed below. 
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4.5 The characteristic equation 

To obtain the solution of Maxwell's equations C lkmn k k k m A n = with in 
the approximation of geometrical optics, we project these equations onto the 



directions X. 



7 (s[T- S (B) b) ) H)(k,A) (B) +A ib) (H (A) -H iB) )} 



J (A) 

Here we introduced the following abbreviations: 

(k, A) (B) = g m ^hA m , = g^kiK 

Explicitly we obtain 

(3) 



-S$<% F = 0. (105) 



(k,A) {B) = TfM 



M (1) M(2)M (3) ^ e {c) k {c) A {c) 



7 



(c)=(l) 



M, 



(3) 



(c)=(l) 



_ H (B) = ± ~ 



/M (1) M (2) M (3) 



(3) 

E 

(c)=(l) 



ill, 



(c) 



The equations l|105(l may then be written as 



siam J ass 



M, 



(a) 



E 



(3) ^ 2 



(c)=(l) 
(3) 



M, 



(c) 



kb) E 



£(c)k(c)A( c 



(c)=(l) 



Mi 



(106) 

(107) 

(108) 
(109) 



+ J 4 (b) i/^ I =0, (110) 



or, in components, 



£ (2) fc (2) g (3) fc (3) A ^ _ £(2)fc(l)fc(2) 



^M (2) M (3 ) M (2) M (3) y 



M (2) 

£(3)fe(l)fe(3) 



(2) 



A/, 



(3) 



A( 3 )=0, (111) 



(i) fc (i) fc (2) A , / , £ (i) fc a) , £ (3) fc ' 3) ^ , 



~M (1) ' \^M (1) M (3 ) ' M (1) ' M (3) 

£(3)fc(2)fe(3) 



M, 



(3) 



A (3) = 0, (112) 
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£(l)fc(l)fc(3) , £(2)^(2)^(3) 



l d) n- — A (2) 



M (1) ^ M (2) 

ffW £ (1) fc2 i) £(2) fc2 2 



fiM (1) M {2) M (1) M (2) 



A (3 ) = 0. (113) 



In addition, we have to consider Eq. I|75|) , which in terms of the tetrad compo- 
nents reads 

k (l)£(l) A (l) + fc (2)£(2)^(2) + fc(3)£(3)-A(3) = 0. (114) 

Formally, we have four equations to find three amplitude coefficients A/ a y How- 
ever, one can check directly, that one of the set of equations (|lllfl - <|113|) . say 
(II 131) . happens to be a linear combination of two others, Ijlllll and Q112I and of 
l)114[l . Thus, we can find one of the quantities, say k(g\E($\A/s), from l|114(l and 
insert it into (|lll(l and (jl 121 This provides us with a system of two equations 
for and A^y. 



H( E ) £ (2)kf 2) s {1) k 2 {1) +e (3) k 2 



(3) 



A(i 



(i) 



^M(2)M( 3 ) M( 2 ) M( 3 ) 

-(4^"43r) £(2)fc(i)fc(2)A2)=0 ' (115) 



1 1 



M (1) M (3) 



'(2) iw (3). 
ff^fc^ £(2)fef 2) +£(3)fc(3) 



AM (1) M (3) M (1) M (3) 



A (2) =0. (116) 



The solution of the system (|115|l . (|116|l requires the Cramer determinant to be 
zero. The resulting characteristic equation is 



H (A) _ H {B) 



± ^jM {1) M (2) M {3) (e {1) k 2 {1) + £ {2 )k 2 2) + £( 3 )kf 3) 



[—-) [( M (2) + M (3)) £ (D fc2 l) + ( M (D + M (3))£(2)^2) 

+ (M (1) + M (a) ) e(3) 4 ) + f-T-J 



0.(117) 



This is a fourth order equation for the wave vector. It corresponds to the 
extended Fresnel equation in (Equations (D.2.23) and (D.2.44)). In the 
isotropic limit Mi = M 2 = M 3 — we have the double solution 

(H^) 2 = (H^) 2 = {g^hkif = 0, (118) 
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which coincides with 1421) . For n 2 = 1 we recover, of course, the vacuum case 



4.6 The uniaxial case 

Obviously, = is also a solution of l|117(l if only two of the three quantities 
Mi, Mi and M3 are zero. Let's assume here M2 = Ms = 0. For the uniaxial 
configuration 



e (2) — £ (3) — e > £ (1) — e (l-£); M(a) — M i 



where < £ < 1, we have 
M 2 = M 3 = 



2 e(l-0 



Mi = £ 



(119) 



(120) 



For the expression multiplying /ft in the brackets in l|117|l we may write 



(M (2) + M (3) ) e (1) /cf 1} + (Af (1) + M (3) ) £(2) fef 2) 

+ (M (1) + M (2) ) e (3) ^ 3) = £e (fcf a) + fcf 3) ) . (121) 



On the other hand, from (|109fl . 

H {A) _ H (B) 



m v 7 

With the help of the last relation the characteristic equation reduces to 



(122) 



H (A) 



(1 



H (B) 







For the choice 



and with the minus sign equations (|123f) factorizes into 

H (A) H (B) = ( g ™(A) kmk ^ (,/,/-•/,,/,. ) Q 



(123) 



(124) 



(125) 



with solutions g mn ^k m k n — and g^^kikj = 0. Both the symmetric tensor 
fields g mn i A ) and g % ^ B ) are optical metrics in such a case, i.e., there are two 
light cones describing the phenomenon of birefringence. The first optical metric 
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IS 



„mn(A) TT'mjjn \rrn \rn ' ( -\rrn \rn i \rm v n \ /ioc\ 

9 =U U - -pX w X (1) - n 2( 1 _ Q \ X (2) X {2) + X (3) X (3) ) ■ ( 126 ) 

The solution g mn ^k m k n — corresponds to a dispersion relation 

where u> = k m U m . This characterizes the so-called extraordinary wave. With 
the refractive index three- vector ri( a ) , defined by 

k( a ) = W7l( ), (128) 



it can be described as a Fresnel ellipsoid 

"(1) n (2) . 4) _ 1 (UQ) 

n 2 + n 2 (l-0 n 2 (l-£) A " UZyj 

In the expression (|104l) for g mn ( B ) the anisotropic contributions exactly cancel 
for 7 = £ _1 and we recover 

mi»(B) = fjmy-n + _L A ™« ( 130 ) 
n 2 



which coincides with the optical metric g* ,nn in i|23|) . Here, the dispersion rela- 
tion is 

^ 2 = ^( fc (l)+ fc ( 2 )+ fc (3)) > ( 131 ) 



which describes ordinary wave propagation. 

Consequently we have obtained an explicit representation (|55|l for the ma- 
terial tensor C lhmn with the optical metrics (|126[> and l|13Ufl . with p, from (|120|l 
and with 7 = £ _1 according to 112411 . This explicit representation within the 
context of the geometric background revealed in (|59[l - 17UI) is the main result 
of this paper. While the general form of the decomposition (|55|l remains true 
in the general case, the interpretation of g mn ( A ) and g mn ( B ) as optical metrics 
is valid only for uniaxial symmetry. Of course, the relations (I59|) - (|62ll between 
the optical metrics g mn ( A ) and g mn ( B ) and g mn W an d g mn ( 2 ) hold in this special 
case as well. This means, a representation (145(1 is possible although g mn W and 
gmn(2) ^ nQ £ navc an obvious physical meaning. 
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4.7 The energy of an electromagnetic wave 



Finally, wc briefly discuss some features of the electromagnetic energy-momentum 
tensor in this approach. The Minkowski tensor of the electromagnetic field reads 

— —QijF rnn H F^ ra H j n g (132) 

Taking into account Maxwell's equations (|33|l and the representation l|32f) of the 
field strength tensor (notice that for the square of the field strength the complex 
conjugate has to be taken), one realizes that the scalar F mn H mn vanishes and 
the energy-momentum tensor takes the simple form 

Tij = 2k{Cj Aik m A n . (133) 

In terms of the optical metrics and with Ij75(l and (|54() the energy scalar is equal 
to 

W = TqlPU* = — -g ln{A) AiA n = -uj 2 e ln AiA n . (134) 
M 

With (|76() . (|79|l and l|81|l one checks that this quantity is positive. Incidentally, 
the energy density of the electromagnetic field does not depend on whether it is 
calculated from the Minkowski or from the Abraham tensor [15] . Remarkably, 
it is only the extraordinary optical metric g ln ( A ) from (|126|l which enters the 
expression (|134|) . The ordinary optical metric g ln ( B ^ does not contribute here. 

Generally, the Minkowski tensor is not symmetric. In the isotropic case, 
however, there exists a generalized symmetry (cf. (E.3.50)). With i) 133(1 . 

1(23), 123) and (001 it takes the form 

glTi = ^ kb k a g* kn A k A n = g* b T\ (135) 
A 4 

in the present context. A corresponding relation is not found in the anisotropic 
case. 



5 Discussion 

We have shown that the material tensor of an anisotropic medium with uniaxial 
symmetry can be constructed out of those two optical metrics which describe 
ordinary and extraordinary light propagation. This generalizes the correspond- 
ing representation of the material tensor in isotropic media in terms of Gordon's 
optical metric. Furthermore, we have clarified the internal geometrical structure 
which underlies the decomposition of the material tensor into a combination of 
products of two symmetric second rank tensor fields. The latter quantities can 
be regarded as vectors in an associated two dimensional space where the coeffi- 
cients of the decomposition play the role of a metric in this space. While we be- 
lieve this type of bi-metric representation in connection with the interpretation 
in terms of an internal metric structure to be new, some already known results 
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are reproduced here in a different context as well. Generally, the wave propaga- 
tion in anisotropic media is governed by a quartic equation for the wave vector, 
called (extended) Fresnel equation This equation is known to factorize 

into two quadratic equations in vacuum. A similar feature is true for isotropic 
media. Here, the factorizing implies a double solution for the wave vector as a 
null vector with respect to the optical metric. For an anisotropic medium with 
uniaxial symmetry a factorization is obtained as well. Instead of the isotropic 
double solution there are two solutions in terms of two optical metrics now, one 
of them describing ordinary, the other one extraordinary wave propagation. All 
these features follow as special cases of our approach. Furthermore, we recover 
a (generalized) symmetry property of Minkowski's energy-momentum tensor for 
isotropic media. 
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